Numerical simulations of conversion to Alfven waves 
in solar active regions 
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__, Abstract. We study the coupling of magneto-acoustic waves to Alven waves using 2.5D 

. numerical simulations. In our experiment, a fast magnetoacoustic wave of a given frequency 

an d wavenumber is generated below the surface. The magnetic field in the domain is assumed 
homogeneous and inclined. The efficiency of the conversion to Alfven waves near the layer 
i-C ■ of equal acoustic and Alfven speeds is measured calculating their energy flux. The particular 

amplitude and phase relations between the oscillations of magnetic field and velocity help us to 
Q . demonstrate that the waves produced after the transformation and reaching upper atmosphere 

• axe indeed Alfven waves. We find that the conversion from fast magneto-acoustic waves to Alfven 

ryj \ waves is particularly important for the inclination 6 and azimuth <j> angles of the magnetic field 

■ between 55 and 65 degrees, with the maximum shifted to larger inclinations for lower frequency 

waves. The maximum Alfven flux transmitted to the upper atmosphere is about 2-3 times lower 
I , than the corresponding acoustic flux. 
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\f} , Conversion from fast-mode high-/3 magneto-acoustic waves (analog of p modes) to slow- 

mode waves in solar active regions is relatively well studied both from analytical theories and 
numerical simulations (e.g., [H [21 13 HI 16] ) , see [7] for a review. In a two-dimensional situation, 
the transformation from fast to slow magnetoacoustic modes is demonstrated to be particularly 
strong for a narrow range of the magnetic field inclinations around 20-30 degrees to the vertical. 
However, no generalized picture exists so far for conversion from magneto-acoustic to Alfven 
waves in a three-dimensional situation. Studies of this conversion were initiated by Cally 
^ & Goossens jS], who found that the conversion is most efficient for preferred magnetic field 

inclinations between 30 and 40 degrees, and azimuth angles between 60 and 80 degrees, and 
that Alfvenic fluxes transmitted to the upper atmosphere can exceed acoustic fluxes in some 
cases. Newington & Cally [9] studied the conversion properties of low-frequency gravity waves, 
showing that large magnetic field inclinations can help transmitting an important amount of the 
Alfvenic energy flux to the upper atmosphere. 

Motivated by these recent studies, here we attack the problem by means of 2.5D numerical 
simulations. The purpose of our study is to calculate the efficiency of the conversion from fast- 
mode high-/3 magneto-acoustic waves to Alfven and slow waves in the upper atmosphere for 
various frequencies and wavenumbers as a function of the field orientation. We limit our study 
to a plane parallel atmosphere permeated by a constant inclined magnetic field, to perform a 
meaningful comparison with the work of Cally & Goossens [8J- Numerical simulation will allow 
generalization to more realistic models in our future work. 

We numerically solve the non-linear equations of ideal MHD assuming all vectors in three 
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Figure 1. Time- height variations of the three projected velocity components corresponding to 
e perp (Alfven wave, left), e tran (fast wave, middle) and ei ong (slow wave, right) for v = 5 mHz in 
a simulation with B inclined by 9 = 30° and <j> = 80°. The solid line marks the position cs = va, 
and the dashed line marks the cut-off layer v = u c cos9. The colour scaling is the same in all 
panels. The amplitudes are scaled with ^Jpqva (first two panels) ^/pocs (last panel). 

spatial directions and all derivatives in two directions (i.e. 2.5D approximation, see [IP} lll|). 
though perturbations are kept small to approximate the linear regime. An acoustic wave of a 
given frequency and wave number is generated at —5 Mm below the solar surface in a standard 
model atmosphere permeated by a uniform inclined magnetic field. The top boundary of the 
simulation box is 1 Mm above the surface, and 0.8 Mm above the layer where the acoustic speed, 
cs, and the Alfven speed, va, are equal. We consider frequencies v = 3 and 5 mHz and wave 
numbers kx = 1-37 Mm -1 and ky = 0. The simulation grid covers field inclinations 6 from 
0° to 80° and field azimuths 4> from 0° to 160°. The field strength is kept at B = 500 G. To 
separate the Alfven mode from the fast and slow magneto-acoustic modes in the magnetically 
dominated atmosphere we use velocity projections onto three characteristic directions: 

eiong = [cos (ft sin 6, sin eft sin 9, cos 9}; 

e P erp = [— cos 4>sm 2 9 sin (j), 1 — sin 2 9 sin 2 eft, — cos 9 sin 6* sin eft]; (1) 

Gtrans = [ — COS 9, 0, COS (ft SU1 9] . 

To measure the efficiency of conversion to Alfven waves near and above the cs = va equipartition 
layer, we calculate acoustic and magnetic energy fluxes, averaged over time: 

F ac = (pivi); F mag = (Bi x (vi x B )//x )- (2) 

Figure [T] shows an example of the projected velocities in our calculations as a function of space 
and time. In this representation the larger inclination of the ridges mean lower propagation 
speeds and vice versa. Note, that by projecting the velocities, we are able to separate the modes 
only in the magnetically dominated atmosphere, i.e. above the solid line in Fig. [TJ The figure 
shows how the incident fast mode wave propagates to the equipartition layer and then splits into 
several components. The Alfven wave is produced by mode conversion above 0.2 Mm (left panel) 
and propagates upwards with the (rapid) Alfven speed, confirmed by almost vertical inclination 
of the ridges. Conversely, the essentially magnetic fast-mode low-/? wave produced in the upper 
atmosphere (middle panel) is reflected, and its velocity variations in the upper layers vanish 
with height. The (acoustic) slow- mode low-/3 wave escapes to the upper atmosphere tunnelling 
over the cut-off layer due to the field inclination of 9 = 30°. The amplitudes of the velocity 
variations of the Alfven wave are comparable to those of the slow wave. 

To confirm the Alfven nature of the transformed waves, as revealed by the projection 
calculations, we checked the amplitude and phase relations for all three modes reaching the 
upper atmosphere. For the Alfven mode the magnetic field B\ and velocity variations V\ should 
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Figure 2. Left panel: Logio of the ratio B\ to Vi/^/JIopo for projected velocities and magnetic 
field variations, averaged over all eft, as a function of 6. Black line: fast mode (e tra n projection); 
red line: Alfven mode (e perp ); blue line: slow mode (ei ong ). Right panel: phase shift between 
the projected variations of V\ and B\, as a function of 9 for selected eft. Red lines: Alfven mode; 
black lines: fast mode. 

be in equipartition (i.e. B\ = Vi/y/flopo), and both magnitudes should oscillate in phase (see 
Priest [I2])- Figure [2] presents the calculations of the amplitude ratio -Bi^//UoPo/^i an d temporal 
phase shift between B\ and Vi, where both velocity and magnetic field variations are projected 
in the corresponding characteristic direction for each mode (Eq. [1]). This calculation confirms 
that, indeed, for all magnetic field orientations 6 and (f>, the amplitude ratio for the Alfven 
mode (e perp projection) is around one (left panel). This is clearly not the case for the slow and 
fast modes. For the fast mode, the amplitude ratio is two orders of magnitude larger, and for 
the slow mode, it is two orders of magnitude lower than one. For the Alfven mode the phase 
shifts group around zero for all <j), unlike the case of the fast mode (right panel). We did not 
calculate the phase shifts for the slow mode as the variations of the magnetic field are negligible. 
Thus, we conclude that the properties of the simulated Alfven mode separated by the projection 
correspond to those expected for a classical Alfven mode. 

An example of the height variations of the acoustic and magnetic fluxes is given in Figure 
[3l The total vertical flux (dotted line) is conserved in the simulations except for the limitations 
caused by the finite grid resolution not resolving slow small-wavelength waves in the deep layers 
(see Fig. Q]). Both acoustic and magnetic fluxes show strongest variations near the conversion 
layer and become constant above it between 0.5 and 1 Mm height. The fluxes reaching the upper 
atmosphere depend crucially on the orientation of the field. In this example, the acoustic flux 
decreases with 6 whilst the magnetic flux increases with 6 and becomes larger than the acoustic 
fluxes for 9 = 60°. As the fast wave is already reflected in the upper atmosphere (see Fig. [[]), 
the magnetic flux at these heights is due to the propagating Alfven wave. 

Finally, Figure gives the time averages of the vertical magnetic and acoustic fluxes at the 
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Figure 3. Examples of the height dependence of the magnetic (solid line) and acoustic (dashed 
line) vertical fluxes, defined by Eq. [21 for v = 5 mHz and several and 0. Solid vertical line 
marks the position cs = va, dashed vertical line marks the cut-off layer v = u c cos9. 
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Figure 4. Vertical fluxes measured at the top of the atmosphere at 1 Mm for waves with v = 5 
mHz (left panels) and 3 mHz (right panels). Upper panels give magnetic fluxes and lower panels 
give acoustic fluxes. 

top of the atmosphere as a function of the field orientation. As proven above, the magnetic 
flux at 1 Mm corresponds to the Alfven mode. At v = 5 mHz, the maximum of the magnetic 
flux corresponds to 9 = 50° and (ft = 65°. This maximum is shifted to larger inclinations 
6 = 65° for waves with v = 3 mHz. The presence of the sharp maximum of the Alfvenic flux 
transmission agrees well with the conclusions made previously by Cally & Goossens [8], though 
the exact position of the maximum is shifted to somewhat larger inclinations. The maximum 
of the transmitted acoustic flux corresponds to inclinations 8 ~ 30° for v = 5 mHz waves, and 
to ~ 55° for v = 3 mHz waves, again, in agreement with previous calculations |8j. The 
absolute value of the fluxes is about 30 times lower for 3 mHz compared to 5 mHz. At some 
angles the Afven magnetic flux transmitted to the upper atmosphere is larger than the acoustic 
flux. However, at angles corresponding to the maximum of the transmission, the Alfven flux is 
2-3 times lower than the corresponding acoustic flux. 

It is important to realize that quantitatively simulating mode transformation numerically is 
a challenge, as any numerical inaccuracies are amplified in such second-order quantities as wave 
energy fluxes. The tests presented in this paper prove the robustness of our numerical procedure 
and offer an effective way to separate the Alfven from magneto-acoustic modes in numerical 
simulations. This will allow us in future to study the coupling between magneto- acoustic and 
Alfven waves in more realistic situations resembling complex solar magnetic structures. 
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